arXiv:astro-ph/0503075vl 3 Mar 2005 


Statefinder diagnostic for coupled quintessence 


Xin Zhang 

Institute of High Energy Physics, Chinese Academy of Sciences 
P.O.Box 918(4), 100039 Beijing, China 

zhangxin@mail.ihep.ac.cn 


Abstract 


The problem of the cosmic coincidence is a longstanding puzzle. This 
conundrum may be solved by introducing a coupling between the two dark 
sectors. In this Letter, we study two cases of the coupled quintessence sce¬ 
nario. (a) Assume that the mass of dark matter particles depends exponen¬ 
tially on the scalar field associated to dark energy and meanwhile the scalar 
field evolves in an exponential potential; ( b ) Assume that the mass of dark 
matter particles depends on a power law function of the scalar field and mean¬ 
while the scalar field evolves in a power law potential. Since the dynamics of 
this system is dominated by an attractor solution, the mass of dark matter 
particles is forced to change with time as to ensure that the ratio between 
the energy densities of dark matter and dark energy becomes a constant at 
late times, and one thus solve the cosmic coincidence problem naturally. We 
perform a statefinder diagnostic to both cases of this coupled quintessence 
scenario. It is shown that the evolving trajectory of this scenario in the s — r 
diagram is quite different from those of other dark energy models. 
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There are more and more evidences [1-3] support that the present universe is dominated 
by dark sectors. Combined analysis of cosmological observations, esp. the Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAP) satellite experiment [2], shows that dark energy (DE) 
occupies about 73% of the energy of our universe, and dark matter (DM) about 23%. The 
usual baryon matter which can be described by our known particle theory occupies only 
about 4% of the total energy of the universe. The accelerated expansion of the present 
universe is attributed to that DE is an exotic component with negative pressure, such as 
the cosmological constant [4] or a scalar field with a proper potential (i.e. the so-called 
quintessence) [5]. The cosmological constant A (or vacuum energy) has the equation of 
state w = — 1. The cosmological model that consists of a mixture of vacuum energy and 
cold dark matter (CDM) is called LCDM (or ACDM). While the so-called QCDM cosmol¬ 
ogy is based upon a mixture of CDM and quintessence field. The energy density and the 
negative pressure are provided by the quintessence scalar field 0 slowly evolving down its 
potential V(0). The equation of state of the quintessence —1 < w < —1/3 is guaranteed 
by the slow evolution. However, as is well known, there are two difficulties arise from all 
of these scenarios, namely, the ’fine-tuning’ problem and the ’cosmic coincidence’ problem. 
The cosmic coincidence problem [6] states: Since the energy densities of DE and DM scale 
so differently during the expansion of the universe, why are they nearly equal today? To get 
this coincidence, it appears that their ratio must be set to a specific, infinitesimal value in 
the very early universe. 

A possible solution to this cosmic coincidence problem may be provided by introducing 
a coupling between quintessence DE and CDM. This coupling is often described by the 
varible-mass particle (VAMP) scenario [7]. The VAMP scenario assumes that the CDM 
particles interact with the scalar DE field resulting in a time-dependent mass, i.e. the mass 
of the CDM particles evolves according to some function of the scalar field 0. In this Letter 
we study two cases of this coupled quintessence scenario: (a) The quintessence scalar field 0 
evolves in an exponential potential and the DM particle mass also depends exponentially on 
0; ( b ) The quintessence scalar field 0 evolves in a power law potential and the DM particle 
mass also depends on a power law function of 0. In both cases, the late time behavior of the 
cosmological equations will give accelerated expansion and, a constant ratio between DM 
energy density p x and DE energy density [8,9]. This behavior relics on the existence of 
an attractor solution, which makes the effective equation of state of DE mimic the effective 
equation of state of DM at late times so that the late time cosmology insensitive to the 
initial conditions for DE and DM. Therefore, the scenario containing coupled quintessence 
with VAMPs solves the cosmic coincidence problem in this sense. 

In this Letter, we will first show the solution to the problem of cosmic coincidence given 
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by this scenario and, then we perform a statefinder diagnostic for both cases of this conpled 
quintessence model. The statefinder parameters introduced by Sahni et al. [10] are proven 
to be useful tools to characterize and differentiate between various DE models. We show in 
this Letter that the evolving trajectory of this scenario in the s — r diagram is quite different 
from those of other DE models. 

Consider, now, the interacting DE model in which we postulate that the DM component 
X interacts with the DE field (j) through the interaction term Q according to 

Px + 3Hp x = -Q , (1) 


p<t> + 3Hp<f,( 1 + w#) = Q , 


( 2 ) 


where p x and p$ are energy densities of DM and DE, respectively, dot denotes a derivative 
with respect to time t, H = a/a represents the Hubble parameter, in which a(t) is the scale 
factor of the universe, and 


Pi = \<t > 2 

P<t> + V(<j>) 


is the usual parameter of equation of state for the homogeneous scalar field <f> associated to 
DE, where V (0) is some potential of the quintessence field. For convenience we can define 
the effective equations of state for DM and DE through the parameters 
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( 4 ) 


If the effective equation of state parameters and evolve to be an equal constant 
at late times, the field system is then proven to have a stable attractor solution. In what 
follows we will discuss two cases of this coupled quintessence scenario — the exponential 
case and the power law case. 

Exponential case 

Assume that the DM particle y with mass M depending exponentially on the DE field 


M x {4>) = M,e~ x<p , (5) 

where <f> is expressed in units of the reduced Planck mass M p ( M p = 1/vo7tG = 2.436 x 
10 18 GeV), and A is a positive constant. The scalar field has an exponential potential 

V(<f>) = Ke^ , (6) 


3 



where p is a positive constant. As a consequence, the interaction term Q in this case can be 
given 


Q X<j>p x • 


(7) 


The equation of motion is then 


1 

3 


°X + Pb + Prad + V) 

I .rt i 




Ap x - Q V , 


( 8 ) 


where pb and p ra( j are energy densities of baryons and radiation, respectively, and we have 
assumed a spatially flat universe. Primes denote derivatives with respect to u — ln(a/ao) = 
— ln(l + z), in which z is the red-shift, and «o represents the current scale factor. Since we 
are interested in the late-time behavior, we can assume pb, p ra d p x , p<y In this limit it is 
easy to see that there is a solution 


0 = 0 0 - 



(9) 


such that 




~ 1 — = 


3 + A(A + p) 
(A + p) 2 


( 10 ) 


and 


wf = wf = W 


X 

A + r/ 


( 11 ) 


This is an attractor in held space for rj > (—A + VA 2 + 12)/2. When the attractor is reached, 
the energy densities of DM and DE will evolve at a constant ratio depending only on A and 
p, thus solving the cosmic coincidence problem. 

It is shown by (11) that W is negative and may lead, if W < —1/3, to an accelerated 
expansion of the universe. To understand how it is possible to get both acceleration and 
constant ratio between DM and DE one may look at the scaling behavior of the energy 
densities on the attractor (9), 


—\<t>—3u r)tj> — 3 ( 1 +W)u 

p x ~ e ~ p<j> ~ e ~ e v ’ . 


( 12 ) 


The scaling behavior of DM deviates from the usual scaling way e~ 3M due to the ^-dependence 
of the DM mass. The interaction between DM and DE forces their effective equation of 
state parameters to become an equal negative constant W, and thus solving the coincidence 
problem and at the same time resulting in an accelerated expansion. 
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FIG. 1. A typical solution to the exponential case. The evolution of the density parameters 
for different components and the effective equation of state parameters for DE and DM. The 
corresponding model parameters are: rj = 2, A = 3, M* = 0.4po/n x o and E* = O.lpo. 


The time evolution of the density parameters for different components (including also pb 
and p r ad) and the effective equation of state parameters for DE and DM for a typical solution 
is plotted in Fig.l. Notice that the attractor solution is going to be reached currently in this 
example. 

Power law case 

In this case we assume that the DM particle x with mass M depending on a power law 
function of the DE field </>, 


M x {4>) = , 

and the scalar field has a power law potential 

V(<t>) = , 


where a,/3 > 0. The interaction term Q in this case is then expressed as 


Q 



The equation of motion can be given 


1 (Px + Pb + Prad + V) It, 
3 l-0' 2 /6 


+ 2 + A> + 2 Prad + 2 V)(f>' 
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-Px 


Jlv 


(13) 


(14) 


(15) 


(16) 
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and it can be proven that there is a stable attractor solution in the held space [9,11] 


0 = 



(17) 


such that 


and 


&,<!> ~ 1 — O x = 


a 

a + f3 1 


(18) 


w^ = W 


a 

ex f3 


(19) 


When the attractor is reached, the energy densities of DM and DE will evolve at a constant 
ratio depending only on a and (3, thus solving the cosmic coincidence problem. We see that 
in this case W is also a negative constant and can thus lead to an accelerated expansion of 
the universe. The scaling behavior of the energy densities on the attractor (17) is exhibited 
as 


„ i—a—'ivu ±0 —3(1 +W)u 

p x ~ <p e ~ Pfj, ~ <p ^ ~ e y ’ 


( 20 ) 



u= -ln(1+z) 

FIG. 2. A typical solution to the power law case. The evolution of the density parameters 
for different components and the effective equation of state parameters for DE and DM. The 
corresponding model parameters are: a = If, j3 = 4, M* = 230po/n x o and K = O.lpo- 
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In Fig.2 we plot a typical solution, including also pb and p ra d- Notice that the attractor 
solution is going to be reached today in this example. 

In what follows we will perform a statehnder diagnostic to this coupled quintessence 
scenario. Since more and more DE models are constructed for interpreting or describing 
the cosmic acceleration, there exists the problem of discriminating between the various con¬ 
tenders. In order to be able to differentiate between those competing cosmological scenarios 
involving DE, a sensitive and robust diagnostic for DE models is a must. For this purpose 
a diagnostic proposal that makes use of parameter pair {r, s}, the so-called ” statehnder”, 
was introduced by Sahni et al. [10]. The statehnder probes the expansion dynamics of the 
universe through higher derivatives of the expansion factor a and is a natural companion to 
the deceleration parameter which depends upon a. The statehnder pair {r, s} is defined as 
follows 

a _ r — 1 

aPP' S = 3(q- 1/2) ' 

The statehnder is a ’geometrical’ diagnostic in the sense that it depends upon the expansion 
factor and hence upon the metric describing space-time. 

Trajectories in the s — r plane corresponding to different cosmological models exhibit 
qualitatively different behaviors. The spatially hat LCDM scenario corresponds to a hxed 
point in the diagram 



!»■ r} 


{ 0 , 1 } . 


( 22 ) 


I LCDM 

Departure of a given DE model from this hxed point provides a good way of establishing the 
’distance’ of this model from LCDM [12]. As demonstrated in [10,12-14] the statehnder can 
successfully differentiate between a wide variety of DE models including the cosmological 
constant, quintessence, the Chaplygin gas, braneworld models and interacting DE models. 
The interacting DE model analyzed in Ref. [14] cannot solve but only alleviate the cosmic 
coincidence problem. We in this Letter will perform a diagnostic for the coupled quintessence 
scenario which can provide a natural solution to the coincidence problem and show explicitly 
the difference between this scenario and other DE models. 

The statehnder parameters can be expressed in terms of the total energy density p and 
the total pressure p in the universe: 


9 (p + p)p 

r — 14-, 

2 P p 


(p + p) i> 
p p ’ 


(23) 


Since the total energy of the universe is conserved, we have p = —3 H(p + p). Then making 
use of p 0 = —3H (1 + w^ e) )p</. and Prad = -4iLp rad , we can get 
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(24) 


p/H = p' = K, - 3^(1 + . 


Hence, the statefinder parameters for the coupled quintessence scenario can be obtained 


r = 1 - - 3ru 0 (l + wJ e) )]H 0 + 2D rad , 


(25) 


-3 [w^ - 3w 0 (l + w^ e) )]H 0 + 4H rad 


s = 


T 3 Q rad 

The deceleration parameter is also given 


(26) 


Q 2 ^rad) • 


(27) 



FIG. 3. The s — r diagram of the exponential case: evolution trajectories of r(s) in the variable 
interval u £ [—2,2], Selected curves r(s) for A = 3, 2, and I, respectively. Dots locate the current 
values of the statefinder pair {s,r}. 


We first apply a statefinder analysis on the exponential case. In Fig.3, we show the time 
evolution of the statefinder pair {r, s}. The plot is for variable interval u £ [—2, 2], and the 
selected evolution trajectories of r(s) correspond to rj — 2 and A = 3, 2 and 1, respectively, 
and the other model parameters are taken to be the same values as those used in Fig.l. We 
see clearly that the distant from this model to LCDM scenario is somewhat far. It is of 
interest to find that the trajectory of r(s) will form swirl before reaches the attractor, which 
is quite different from other DE models (see [10,12-14]). It is demonstrated again that the 
statefinder can successfully characterize and differentiate between various DE models. As 
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complementarity for the diagnostic, we also plot the evolution trajectories of statehnder pair 
{r, q} in Fig.4. We see that the cosmic acceleration is ensured by sufficient strong coupling 
A. 



FIG. 4. The q — r diagram of the exponential case: evolution trajectories of r(q) in the variable 
interval u E [—2,2], Selected curves r(q) for A = 3, 2, and 1, respectively. Dots locate the current 
values of the statehnder pair {g, r}. 



FIG. 5. The s — r diagram of the power law case: evolution trajectories of r(s) for the variable 
interval u E [—2,2], Selected curves r(s) for a = 11, 13, and 15, respectively. Dots locate the 
current values of the statehnder pair {s,r}. 
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Next we apply a statefinder diagnostic to the power law case. In Fig.5, we show the time 
evolution of the statefinder pair {r, s}. The plot is also for variable interval u G [—2,2], 
and the selected evolution trajectories of r(s) correspond to (3 = 4 and a = 11, 13 and 
15, respectively, and the other model parameters are as the same as those used in Fig.2. 
It can be seen that the trajectories of this case will pass through LCDM fixed point. And 
the swirls in this case are more evident than those of exponential case. We also plot the 
evolution trajectories of statefinder pair {r, q} in Fig.6. 
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FIG. 6. The q — r diagram of the power law case: evolution trajectories of r(q) for the variable 
interval u 6 [—2,2]. Selected curves r(q) for a = 11, 13, and 15, respectively. Dots locate the 
current values of the statefinder pair {q, r}. 


In summary, we study in this Letter the statefinder of the coupled quintessence scenario. 
We analyze two cases of this scenario — the exponential case and the power law case. It 
is shown that both cases of this scenario have attractor behaviors and can thus provide a 
natural solution to the cosmic coincidence problem. Then we perform a statefinder diagnostic 
to both cases of this coupled quintessence scenario. It is shown that the evolving trajectory 
of this scenario in the s — r plane is quite different from those of other DE models. We hope 
that the future high precision observations (e.g. SNAP) will be capable of determining these 
statefinder parameters and consequently shed light on the nature of DE. 
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